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Abstract 



The possibility to study elastic and diffractive physics in pp collisions at LHC is dis- 
cussed. For this purpose we have considered detectors close to the beam in conjunction 
with the magnetic elements of the accelerator to provide a high precision spectrometer for 
very forward final state protons. The geometrical acceptance is given and the momentum 
resolution is calculated for different spatial resolution detectors. 



1 Introduction 



Elastic and diffractive scattering processes are generally interpreted in terms of Pomeron 
(IP) exchange . They are characterized by a final state which contains particles emitted 
at very small angles with respect to the incoming beam and carry a substantial fraction 
of the available energy. At collider experiments these particles escape undetected down 
the beam pipe. The Roman pot technique @ allows to position detectors directly in the 
vacuum pipe at a very small distance from the beam. 

Since 1960 elastic and diffractive scattering processes have been studied in hadron- 
hadron interactions at CERN || and elsewhere. Diffractive events in deep inelastic ep 
scattering have been also measured at HERA H. Within the ZEUS experiment, the 
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Leading Proton Spectromenter (LPS) |5j has been built to measure the protons scattered 
at very small angle with respect to the incoming beam direction. 

In this paper we investigate the possibility to study elastic and diffractive events in pp 
collisions at LHC energy (y/s = 14 TeV). In particular we study single diffractive events 
characterized by one of the protons emerging from the interaction isolated in rapidity 
and with its momentum very close to that of the beam while the other diffracts into a 
state of mass Mx- The Feynman graphs of the interaction are shown schematically in 
Fig.[TJ. Preliminary studies on this process are presented in 0. A full acceptance detector 
could be built in the last free LHC intersection region (IP4) where one could have low 
luminosities, of the order of 10 30 cm~ 2 s _1 . An experimental apparatus for the complete 
forward region has been discussed in 0, and JTJ . 



In the following we present the simulation of spectrometers in the very forward region 
along one of the two outgoing proton beam arms at a sufficient distance from the inter- 
action point in order to accept final state particles produced at very small angle with 
respect to the incoming proton beam direction. The present paper is an extended and 
improved version of ||. 



2 Simulation of elastic and diffractive pp events 



The process under study is the reaction pp — > pX (Fig.[IJ) in which one of the protons 
remains intact and loses a small fraction of its momentum 1 — xl, with xl = ^| , where 
p and p' are the incoming and outgoing proton momenta, respectively. This quantity is 
related to the mass Mx of the hadronic system produced in the interaction by the relation: 



£ = i - x L 



Ml 



s 

where s is the square of total centre of mass energy. Elastic scattering is characterized by 
xl = 1 while diffractive processes dominate for xl > 0.9. The transverse momentum 
of the scattered proton is generally small: pt < 1 GeV/c. 

In order to select the possible regions along the proton beam line where these processes 
can be studied, we need to examine in detail the properties of the proton beam and the 
geometrical restrictions imposed by the vacuum pipe. 

Two possible ways to study the beam profile as it proceeds through the magnetic 
elements are summarized in the following. The first method makes use of the beta function 
(3 and the phase advance A// parameters extracted directly from the LHC MAD tables 
|T0|j . The transverse distance in the bending plane of a particle from the nominal beam, 



x, and its angle 9 X with respect to the nominal beam direction can be expressed as: 
x{l) 



\ &(0) 



cosAfix(0) + J p x (l) (3 X (0) smAfi 6 X (0) + £ D x , (1) 
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and 



9 x (l) = -, ^=(sinAfi+a x (l)cosAfi)x(Q) + 



ir{ °' {cosAfi-a x {l)sinAn)6 x {0)+£D' x , 



(2) 

where 

a,(0 = "0.5 ^ (3) 

and D is the dispersion. In the case of elastic scattering £ is the natural beam momentum 
spread of 10~ 4 . 

The position and angle of the particle at the interaction point are indicated as x(0), 9 X (0). 
Here I is the coordinate along the beam path. Similar equations can be written for the y 



plane. Results using this method are given in || and p~3 [ . 

The second way to look at the beam profile is to use the transport matrices of the 
beam line: 




H21 H22 




(4) 



where the matrix elements Hy are functions of xl only. 

Fig.||| shows a possible layout for an insertion in intersection region IP4 at LHC which is 
consistent with the constraints imposed by the experimental magnets and RF cavities M 



and [bj] . The experimental equipment for the detection of very forward scattered protons 
may be located in Roman pots inside the vacuum beam pipe. Tables list the main 
features of the magnetic elements of one of the outgoing arms, from the interaction point 
up to 450 m. 

Fig.|3| shows the 10a profiles of the beam in the horizontal and vertical planes. They 
have been calculated using both methods discussed earlier. The orbits agree all along the 
beam line, showing that the tracking methods are equivalent. 

All the results presented in this paper have been obtained with a program which sim- 
ulates particle trajectories using the beam optics transport equation (^). This program 
reflects in its main features the TRANSPORT beam optics program |lTJ and is similar to 



the simulation program used to design the Leading Proton Spectrometer at ZEUS 0, |T3] . 
The beta value at the intersection point considered in the calculation is 1000 m. 

A short description of the main characteristics of the simulation program is given in the 
following. Given the beam momentum, a track is followed from the starting point until it 
is lost at an obstruction. The program returns the hit pattern at selected positions where 
detectors may be placed. The elements along the beam line which have been simulated 
are: 

• Drifts - in which no particle is lost and which have the same aperture as that of the 
previous element. 

• Quadrupoles - the focussing effect is calculated as a function of xl- Negative focal 
length means horizontal focussing. Negative xl means particle with negative charge. 
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• Dipoles - the bending power is calculated with respect to that for a track at xl = 
1, so that no bending occurs for particles at Xl = 1. During tracking inside a 
magnetic element the distance of the particle from the nominal beam may reach a 
local maximum and a test is done to see if this maximum is outside the beam pipe. 
The same test is performed at the entrance and at the exit of each quadrupole. 

• Detectors - at present these are considered to be virtual planes covering the whole 
beam pipe aperture except for a rectangular hole with sides of length equal to I0a x 
and 10a y , where a x and a y are the horizontal and vertical standard deviation of the 
beam distribution (taken to be a Gaussian), respectively. 

For the following study particles were generated in the interaction vertex diamond of 
dimensions (Ax, Ay, Az) =(0.6mm, 0.6mm, 0.0mm), according to a flat distribution in 
log(l — Xl) within the limits 10~ 4 < (1 — xl) < 10 _1 and a flat distribution in log{px) 
within the limits 10~ 4 < pr < 1 MeV/c. The z axis points in the proton beam direction. 
The Xl distribution has been chosen because it reflects the observed single diffractive 
cross section, the px distribution in order to enhance the generation of particles at very 
low pt (and thus at very low \t\) and therefore approach the optical point (t = 0) in the 
study of the elastic cross section at (1 — xl) = 0. 



3 The spectrometers 

In order to maximize the geometric acceptance all the detectors are placed in the region 
beyond 200 m from the IP where the beam profile has the smallest transverse size (see 
Fig-D- 

In this way we select three regions where the beam elements can be used to momentum 
analyse very forward outgoing protons. We call them spectrometers A, B, C. Each of 
them consists of two detector planes positioned upstream and downstream of a dipole ele- 
ment. The detector planes are located only where the distance between adjacent magnetic 
elements is > 50 cm. 

Fig.|] shows particle trajectories through one of the outgoing arms from IP4. They are 
calculated considering two values of xl and three values of pr- We choose two extreme 
Pt values of ±2 GeV/c and the central value pt = 0. The beam pipe diameter has been 
taken to be infinite when tracking particles at this stage. Vertical bands are superimposed 
to indicate the positions of the three spectrometers. 

Table |] summarizes the positions and characteristics of the planes of each of the spec- 
trometers. 

In addition to spectrometers A, B, C which cover the very forward region, there is a 
45 m long drift space downstream of the central experimental region between dipoles Dl 
and D2 which allows the detection of the forward going leading particles with low xl down 
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to xl ~ 0. The geometrical acceptance is calculated also for this region in order to see 
which range of xl and pt can be detected outside the vacuum beam pipe. A preliminary 
study of this region is given in 0. 



4 Geometrical acceptance and momentum resolution 

At least three points are needed to measure the momentum. We use the interaction 
vertex and demand the coincidence of any two detector planes (double coincidence events). 
Using the transport matrix (|j) we relate the impact point of the track on the detector to 
the position and angle of the track at the interaction point. At plane i we can write the 
equation 

Xl = H m x{0) + H i;12 9 x {0). (5) 

Assuming the nominal vertex (x(0)=0) we can also relate the position of the track at two 
different planes i, j with the equation 

Xi = Mij^x^xj + C i:j (x L ). (6) 

For a given pair of detectors i, j and for a given value of xl, the coordinates Xi and Xj 
are thus linearly correlated; the parameters of the straight line are functions of xl only, 
thereby allowing the determination of xl from double coincidence events. 

This method is used for momentum reconstruction in the ZEUS LPS [||] and a lon- 
gitudinal momentum resolution of 0.4% has been achieved for particles at xl = 1; the 
resolution on the transverse plane for the ZEUS LPS is less than the beam transverse 
momentum spread of ~ 40 — 100 MeV/c. 

Fig.|5| shows the pattern of hits at the detector planes of the three spectrometers. The 
10(T profile of the beam is often clearly visible. The area covered by the hits on each plane 
is always asymmetric. This asymmetry originates from the bending and focalizing power 
of all the magnetic elements of the proton beam line upstream of the detector planes; this 
must be taken into account in order to define the best geometry of the detector planes. 

Fig.|6| shows the correlation between horizontal and vertical coordinates for the different 
pairs of detector planes which are used to calculate the momentum. The regions of the 
planes which detect very high xl and very small pr tracks correspond, in the figure, to 
the areas of higher event density. This is because of the particular generation we have 
used which is mentioned in the previous section. The straight lines, superimposed to the 
scatter plot shown in Fig.||c, correspond to the different values of xl at which particles 
were generated. Correlations on the horizontal plane provide a good xl resolution because 
the lines which represent two adjacent values of xl are well separated. On the other hand 
on the vertical plane the lines representing the same values of xl would be densely packed 
and are not used for resolution calculation. 
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The geometric acceptance for the double coincidence events is shown in Fig.fT], Fig.|] 
and Fig.|9| as a function of £ and px for the spectrometers A, B and C, respectively. 
The acceptance is limited by the vacuum beam pipe apertures at the highest values of £ 
and the highest values of px- On the other hand particles with very small px and with 
momentum close to the beam momentum escape within the 10<7 detector apertures. 

From Fig.0, Fig.[8] and Fig.[| we can conclude that spectrometer A has 100% acceptance 
in the region 0.03 < £ < 0.1 and 1 < px < 10 3 MeV/c and in the region 10~ 4 < £ < 0.1 
and 10 2 <Pt% 10 3 MeV/c. Spectrometer B has 100% acceptance in the region 10~ 2 < 
i < 0.03 and 1 < p T < 10 3 MeV/c and in the region 10" 4 < £ < 0.03 and 50 < p T < 500 
MeV/c. Finally spectrometer C has 100% acceptance in the region 10~ 3 < £ < 10~ 2 and 
1 % Pt % 10 3 MeV/c. In general very low £ tracks are accepted for 10 2 %, Px H 10 3 
MeV/c, and very low p T tracks are accepted in the region 10~ 3 < £ < 0.1. 

A different region of £ can be studied by a spectrometer placed in the 45 m long 
drift region between dipoles Dl and D2. Two detector planes have been assumed at the 
beginning and at the end of this drift space and their geometrical acceptance is shown in 
Fig. 10 as a function of £ and px- In this case £ has been generated uniformly between 
and 1. The figure shows a wide interval of £ between 0.15 and 0.6 in which the acceptance 
at the beginning of the drift space is 100% for particles inside the pipe. The acceptance 
does not depend on p T . At the end of the drift space only particles with £ < 0.2 stay 
inside the beam pipe and, if we require the coincidence of the two planes, particles are 
accepted only in a very limited region of £. This means that particles with 0.2 < £ <^ 0.6 
can be analysed with detectors (and magnets) placed outside the vacuum pipe. 

The momentum resolution has been calculated for the three spectrometers with 
either 100/im or 10/xm pitch detectors and is shown in Fig.|TT|. Two different sets of 
curves are shown in the figure. The continuous lines show the momentum resolution for 
the three spectrometers A, B, C equipped with 100/im pitch detectors, while the dashed 
lines show the resolutions for 10/im pitch detectors. All curves show an approximately 
constant behaviour in the considered (1 — xl) interval, and in general the spectrometers 
with 10/mi pitch detectors show a resolution which is one order of magnitude better than 
the stations with 100/mi pitch detectors. The resolution improves as the distance from the 
interaction point increases. This effect originates from the increased integrated magnetic 
field acting on the particle. The best resolution is obtained with spectrometer C and has 
the values ^ ~ 10~ 3 % and ^ ~ 10~ 4 % for the lOOtxm and lOum pitch detectors, 
respectively. Spectrometer C accepts however only particles with (1 — x^) < 10~ 2 due to 
the large distance of the detectors from the interaction point. We can conclude that the 
order of magnitude of the momentum resolution is always better of 1 GeV/c for a typical 
momentum of 7 ■ 10 3 GeV/c. 

The momentum resolution calculated using vertical coordinate correlations is very poor 
compared to the previous one. This can be seen also from Fig|| where the vertical 
correlation lines are very densely packed and generate many ambiguities in the momentum 
calculation. 

The transverse momentum resolution — (PT(rec)-p T (gen)) ag f unc tion of pxiQen) 
has been studied for the three spectrometers separately, with either 100/im or 10/im pitch 
detectors. 
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The transverse momentum px has been calculated using the two components p x and 
p y computed using the relation: 



Px(0) = &x(0) X L (hor) Pbearn, (?) 

where 9 X (0) is obtained from the transport matrix eq.(^), and xi{hor) is the momentum 
estimate in the horizontal plane. The p y component is calculated using a similar equation 
in the y plane, but still considering the momentum estimate XL(hor). 

Fig.[12| shows the pt resolution for spectrometer A for two different detector pitches. 
We can observe that the resolution is always better that 3% for px > 200 MeV/c and 
slowly improves as a function of px- No xl dependence is observed both in the case of 
100/im and 10/im pitch detectors. The resolution is a factor of two better for the 10/im 
pitch detectors planes in all px bins considered. The best resolution in spectrometer A 



is of the order of 0.1%. Fig.[l3| shows the pr resolution for the spectromenter B and C. 
Full and open symbols represent the resolution for the 100/im and for the 10/im pitch 
detectors, respectively. In all cases the resolution slowly improves as a function of px and 
takes values which are approximately of the same order of magnitude as those observed 
in spectrometer A. Finally it is interesting to notice that the calculated resolution in the 
whole Pt range considered is of the order of the nominal beam transverse momentum 
spread (< 5 MeV/c). 



5 Conclusions 



Geometrical acceptance and momentum resolution have been calculated for three spec- 
trometers along one of the outgoing proton beam arms of IP4 at LHC in the interval 
10~ 4 ^,0-1 and <pr ^ 10 3 MeV/c. The apparatus has been simulated with 100/im 
or with 10/im pitch detector planes. We have shown that the three spectrometers cover 
regions of xl which partially overlap. They have very different momentum resolutions 
for equal values of xl but their pt resolution do not differ dramatically. In general going 
from 100/im to 10/im pitch an order of magnitude in the momentum resolution is gained, 
while only a factor two is gained for the px resolution. The momentum resolution is in 
the range 10" 4 % < ^ < 5 ■ 10~ 2 %. The p T resolution is in the range 0.1% < ^ < 3% 
for particles with > 200 MeV/c. 
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Name 


Type 


start 


length 


spot considered 


parameter 






[m] 


[m] 


[mml x [mml 

L J L J 


[mradl if bend 












[ml if quad 

L J J- 


LALEPH 


drift 


0. 


3.83 






DUA1 


bend 


3.83 


4.95 


_ 


0.170 


L2 


drift 


8.78 


2.72 






DO 


bend 


11.5 


2.0 




0.342 


L3 


drift 


13.5 


4.5 






Dl 


bend 


18.0 


9.45 


44x44 


1.62 


L4 


drift 


27.45 


0.6 






Dl 


bend 


28.05 


9.45 


44x44 


1.62 


L5 


drift 


37.5 


45.265 






D2 


bend 


82.765 


9.45 


44x44 


-1.62 


L6 


drift 


92.215 


0.6 






D2 


bend 


92.815 


9.45 


44x44 


-1.62 


L7 


drift 


102.265 


2.2 






Ql 


quad 


104.465 


5.5 


35 x 35 


74.822 


L8 


drift 


109.965 


4.5 






Q2 


quad 


114.465 


5.5 


35 x 35 


-52.397 


L9 


drift 


119.965 


1.0 






Q2 


quad 


120.965 


5.5 


35 x 35 


-52.397 


L10 


drift 


126.465 


5.5 






Q3 


quad 


131.965 


5.5 


35 x 35 


39.440 


Lll 


drift 


137.465 


1.8 







Table 1: Positions of the beam elements of one of the two outgoing proton beams with 
respect to the IP4 and the corresponding field values considered in the simulation program. 
The spot considered is the area transverse to the nominal beam available for the passage 
of particles. 
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Name 


Type 


start 


length 


spot considered 


parameter 






[m] 


[m] 


[mm] x [mm] 


[mrad] if bend 












[m] if quad 


L12 


RF cav 


139.265 


33.528 






Q4 


quad 


172.793 


3.0 


28 x 28 


-123.001 


L13 


drift 


175.793 


20.0 






Q5 


quad 


195.793 


3.0 


28 x 28 


85.251 


L14 


drift 


198.793 


1.2 






D3 


bend 


199.993 


9.45 


38 x.38 


-1.62 


L15 


drift 


209.443 


0.6 






D3 


bend 


210.043 


9.45 


38 x.38 


-1.62 


L16 


drift 


219.493 


16.0 






D4 


bend 


235.493 


9.45 


38 x.38 


1.62 


L17 


drift 


244.943 


0.6 






D4 


bend 


245.543 


9.45 


38 x.38 


1.62 


L18 


drift 


254.993 


1.0 






L19 


drift 


255.993 


0.2 






Q6 


quad 


256.193 


3.0 


28 x 28 


-45.045 


L20 


drift 


259.193 


0.2 






Q6 


quad 


259.393 


3.0 


28 x 28 


-45.045 


L21 


drift 


262.393 


2.5 






BDSO 


bend 


264.893 


14.2 




5.1 


L22 


drift 


279.093 


1.46 






BDSO 


bend 


280.553 


14.2 




5.1 


L23 


drift 


294.753 


2.14 






Q7 


quad 


296.893 


3.25 


22 x 18 


34.150 


L24 


drift 


300.143 


0.2 






Q7 


quad 


300.343 


1.5 


22 x 18 


1536.098 



Table 2: Same as Table 1. 
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Name 


Type 


start 


length 


spot considered 


parameter 






[ml 

L J 


[ml 

L J 


[mm] 


[mrad] if bend 












[m] if quad 


L25 


drift 


301.843 


8.91 






BDSO 


bend 


310.753 


14.2 




5.1 


L26 


drift 


324.953 


1.46 






BDSO 


bend 


326.413 


14.2 




5.1 


L27 


drift 


340.613 


2.14 






08 


quad 


342.753 


3.25 


22 x 18 


-34.150 


L28 


drift 


346.003 


0.2 






Q8 


quad 


346.203 


1.5 


22 x 18 


-139.762 


L29 


drift 


347.703 


4.86 






BDSO 


bend 


352.563 


14.2 




5.1 


L30 


drift 


366.763 


1.46 






BDSO 


bend 


368.233 


14.2 




5.1 


L31 


drift 


382.423 


2.14 






Q9 


quad 


384.563 


3.25 


22 x 18 


34.150 


L32 


drift 


387.813 


0.2 






Q9 


quad 


388.013 


1.5 


22 x 18 


175.438 


L33 


drift 


389.513 


4.99 






BDSO 


bend 


394.503 


14.2 




5.1 


L34 


drift 


408.703 


1.46 






BDSO 


bend 


410.200 


14.2 




5.1 


L35 


drift 


424.363 


9.321 






Q10 


quad 


433.684 


1.5 


22 x 18 


-959.233 


L36 


drift 


435.184 


0.2 






Q10 


quad 


435.384 


1.5 


22 x 18 


-73.992 



Table 3: Same as Table 1. 
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Spectrometer A 


Detectors 


Position 


a x of beam 


o y of beam 




[m] 


[mm] 


[mm] 


Al 


175.793 


0.50 


0.27 


A2 


219.493 


0.37 


0.01 


Spectrometer B 


Detectors 


Position 


o^of beam 


o y of beam 




[m] 


[mm] 


[mm] 


Bl 


301.843 


0.17 


0.06 


B2 


340.613 


0.04 


0.04 


Spectrometer C 


Detectors 


Position 


o^of beam 


Oy of beam 




[m] 


[mm] 


[mm] 


CI 


389.513 


0.23 


0.09 


C2 


424.363 


0.13 


0.21 



Table 4: List of detector planes considered in the three spectrometers and la beam 
at the planes. 
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Figure 1: Single pomeron (IP) exchange in elastic (a) and single diffractive (b) pp inter- 
actions. In interactions of type (b) a final hadronic system Mx is produced. 
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IP Dl D2Q1 Q2Q3RFQ4Q5D3 D4 Q6 Q7 Q8 Q9 Q10 




Figure 2: Plan view of a possible insertion in IP4 showing the positions of the three 
spectrometers considered in this paper (from [8]). 
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Figure 3: 10a profiles a) horizontal; b) vertical. The continuous lines show the orbits 
calculated with eqs.(l),(2) (as in [8]) the dashed lines those calculated with eq.(4). 
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Figure 4: Typical horizontal and vertical trajectories of particles through the magnetic 
elements of the beam line at different Xl and p? values. The continuous lines correspond 
to particles with px = 0, the two dashed lines correspond to particles with pt = ±2 
GeV/c. 
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Figure 5: Distribution of the hits at the detector planes of spectrometers A, B and C. A 
box is drawn at each plane representing the dimensions of the 10a profile of the beam. 
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Figure 6: Correlations between track positions in pairs of planes of the spectrometers. 
Continuous, dashed, dot-dashed lines in (c) correspond to tracks with xl — 0.999,X£ = 
0.980,xl = 0.945, respectively. 
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Figure 7: Geometrical acceptance as a function of £ and pr for double coincidence events 
in spectrometer A. 
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Spectrometer B 
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Figure 8: Geometrical acceptance as a function of £ and pr for double coincidence events 
in spectrometer B. 
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Spectrometer C 
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Figure 9: Geometrical acceptance as a function of £ and pr for double coincidence events 
in spectrometer C. 
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Figure 10: Geometric acceptance as a function of £ and pr for two planes placed at the 
beginning (a) and at the end (b) of the 45 m drift space between dipoles Dl and D2. 
Three acceptance ranges are shown in the pictures: ~ 100% (dark grey), 100 — 50%, and 
50 - 0% (light grey). 
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Figure 11: Momentum resolution ^— for the spectrometers A, B and C. Results for 
100/xm (continuous lines) and 10/xm (dashed lines) pitch detectors are shown. 
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Figure 12: Transverse momentum resolution for spectrometer A for particles gener- 
ated in two different Xl bins and for 100/xm (full symbols) or 10/xm (open symbols) pitch 
detectors. 
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Figure 13: Transverse momentum resolution for spectrometers B and C equipped 
with 100/xm (full symbols) or 10/xm (open symbols) pitch detectors. 
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